Proper ferroelastic phase transitions in thin epitaxial films with symmetry-conserving 

and symmetry-breaking misfit strains 
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We study how the ferroelastic domain structure sets in in an epitaxial film of a material with 
second order proper ferroelastic transition. The domain structures considered are similar to either 
a.\l aij a\l ai or c/a/c/a structures in perovskite ferroelectrics. If the "extrinsic" misfit strain, not 
associated with the transition, does not break the symmetry of the high-temperature phase, the 
phase transition in the film occurs at somewhat lower temperature compared to the bulk. The 
loss of stability then occurs with respect to a sinusoidal strain wave, which evolves into the domain 
structure with practically the same geometry and approximately the same period. In the presence of 
the symmetry-breaking component of the misfit strain ( "extrinsic" misfit) the character of the phase 
transition is qualitatively different. In this case it is a topological transition between single-domain 
and multi-domain states, which starts from a low density of the domain walls. 

77.80.Dj, 77.55.-ff, 81.30.Dz 



I. INTRODUCTION 

Symmetry breaking spontaneous strains appear fre- 
quently during phase transitions in various materials, 
and. these materials behave as either proper or improper 
ferroelastics. Ferroelastic phase transitions in thin films 
usually lead to formation of domain structures in order to 
accommodate so-called "misfit stresses" , which appear at 
the transition or may already exist in the high-symmetry 
phase because of the difference in the unit cell parameters 
of a film and a substrate |Q| . Currently there is much in- 
terest in the ferroelastic domain structures in thin films of 
"pure" ferroelastics (see, e.g. [||J^]) and perovskite ferro- 
electrics, which are improper ferroelastics (see references 
in Ref. §). 

The continuous-medium description applies to these 
domain patterns, since usually the thickness of the films 
is by two to three orders of magnitude larger than the 
unit cell parameter. One can also expect that the results 
of this approach will be qualitatively valid even for very 
thin films. It should be noted that a rigorous continu- 
ous medium treatment of ferroelastic domain structure 
is a problem of both non-linear and non-local elasticity 
in inhomogeneous medium (e.g. film on a substrate) and 
presents tremendous mathematical difficulties. Even af- 
ter being reduced to a linear local elasticity problem in a 
homogeneous medium (i.e. for the same elastic constants 
of the film and the substrate) it is usually approached 
with a combination of approximate and numerical meth- 
ods (see, e.g. |§-^). Recently we have treated two simple 
yet generic examples of domain structures, which corre- 
spond to either ai/a2/ai/a2 or c/a/c/a domain patterns 
in perovskite ferroelectrics llG,0. Unlike the previous 



authors, we were able to obtain the analytical expres- 
sions for energies of the domain structures. This allowed 
us to reproduce the available results, and obtain the new 
ones. In the present paper we shall analyze the loss of 
stability of the symmetric phase, which eventually leads 
to formation of the domain pattern. 

We shall start with the question of how the domain 
structures appear at the phase transition |]l^. More 
specifically, we suppose that there occurs a second-order 
proper ferroelastic transition in a free sample, charac- 
terized by a one-component order parameter. We con- 
sider first the case where the substrate induces no misfit 
strain apart from that appearing at the phase transition. 
If this misfit strain can be accommodated by a domain 
structure, the domain pattern inevitably appears at such 
a transition. It proceeds with the loss of stability of 
the symmetric phase with respect to a "strain wave" , 
which then evolves into the domain structure. The pro- 
cess is somewhat different for an improper ferroelastic: 
the domain structure there may not appear immediately 
at the transition, since there a single-domain state is 
metastable. Therefore, the domain structure may de- 
velop there at lower temperatures. However, the struc- 
tures are similar in both proper and improper ferroe- 
lastics not very close to the critical point, since they 
are mainly defined by the symmetry of the spontaneous 
strain. Therefore, theoretical study of proper ferroelas- 
tic transitions in thin epitaxial films is of interest on its 
own merits, and because it reveals the general proper- 
ties of domain structures, even when the strain is not a 
"primary" order parameter. 

We find the form of the "strain waves" , which may 
appear at the second order transition. The problem is 
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non-local but linear, that is why we cannot find the am- 
plitude of the "frozen strain wave" just below the tran- 
sition. However, the form and the period of the "strain 
wave" close to the transition proves to be, in some cases, 
practically the same as far from the transition, where this 
"wave" transforms into the domain structure. Therefore, 
the important features of domain structures can be found 
within a linear theory, and this makes the problem of the 
stability loss even more interesting. 

The misfit strain, not associated with the phase tran- 
sition, if any, changes the situation significantly if it con- 
tains the strain corresponding to the order parameter. 
Then, the phase that has been symmetric in the free sam- 
ple is no longer symmetric in the epitaxial film. However, 
a continuous phase transition may still take place. The 
homogeneous non-symmetric state, which does not "feel" 
the lowering of temperature below the phase transition 
in the free sample, may transform, at yet lower temper- 
ature, into an inhomogeneous multidomain state. This 
transition is similar to commensurate-incommensurate 
transition or the transition between homogeneous and 
vortex states in type H superconductors in magnetic field. 
Such transitions, while being continuous, are not accom- 
panied by a stability loss. The phase diagrams (including 
those in the plane of variables misfit strain-film thickness) 
of transition between single- and multidomain states were 
studied by several authors Some properties of the 

domain structure near this transition we have discussed 
earlier We briefly discuss this topic below in the 

context of the present paper. 

We shall consider the epitaxial film with the plane per- 
pendicular to z— axis. The misfit stress, not associated 
with the phase transition ( "extrinsic" misfit), contains 
two components of the strain tensor: u^,^ = \uyy \ , with 
the unstrained state for the free film being just above 
the phase transition. In Sec. 2 we analyze how the do- 
main structure sets in at a proper ferroelastic transition 
with Uxy the order parameter. We find the period of the 
"strain wave" appearing at the transition and estimate 
the critical thickness of the film where the ferroelectric 
phase exists. The extrinsic misfit is not important in this 
case. In Sec. 3 we treat the case of u^x ~Uzz a-s the order 
parameter for the same film geometry, first for the case 
when the permanent misfit is absent and then when it is 
present. The present results are summarized at the end. 



II. STRAIN Uxy AS THE ORDER PARAMETER 

Consider the substrate as an infinite isotropic medium 
with the shear modulus The elastic moduli of the fer- 
roelastic are supposed to be the same as in the substrate, 
with the exception of the "soft" modulus corresponding 
to the Uxy component of strain. Fig. 1. One can guess 
the important features of the ferroelastic domain struc- 
ture by just considering the symmetric phase. Indeed, 
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FIG. 1. Schematic of the a\/a2/ai/a2 "strain wave" (top 
panel) and the domain structure (middle and bottom panels) 
in thin epitaxial ferroelastic film with the spontaneous strain 
ul-y. The symmetric phase looses its stability at some tem- 
perature Tci below the phase transition temperature in free 
film, Tc, and forms the "strain wave" (top panel). Upon fur- 
ther coohng it transforms into the domain structure with the 
period 2a (bottom panel) in the film of thickness I. 

the domain pattern evolves from a "frozen strain wave" 
appearing at the transition. Given the macroscopic scale 
of the wave, one expects that the corresponding strains 
are mainly the "soft" ones. This immediately provides 
the possible directions of the "wave front" , which corre- 
sponds to the possible orientations of the domain walls. 
In the present case the front is perpendicular either to x- 
or y-axes, since for all other directions the acoustic modes 
are not "soft" . Indeed, they involve not only the soft Uxy- 
component of the strain tensor but also other ( "hard" ) 
components. Therefore, to study the stability loss we 
may consider either the waves of the displacement vector 
component Ux , propagating along the y-axis or swap the 
y— and x— axes. Evidently, each wave is characterized by 
the strain 2uxy = (dux/dy) + {duy/dx) and the rotation 
"^^xy = (dux/dy)- (duy/dx) with \flxy\ = \uxy\ ■ This 
elementary result illustrates the fact that a ferroelastic 
domain structure involves both strains and rotations, as 
discussed at some length in, e.g., Ref. |l|l. The wave we 
consider is confined within the film, therefore, we have 
to take into account the dependence of Ux (or Uy) on z. 

It is obvious from the above that we have to retain in 
the elastic (free) energy the terms dependent on Uxy and 
Uxz (or Uxz) only. Among the gradient terms there are, 
strictly speaking, those containing the gradients of flik 
llisl , but their account does not affect the results and 
we shall omit them together with the gradients of the 
"hard" components of the strain tensor. Therefore, we 
choose the Landau free energy in the form 
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where A = a{T — Tc), with a^D^fx positive constants. 
Recall that higher order terms are not needed for the 
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analysis of loss of stability of the symmetric phase. We 
obtain, therefore, only one non-local equation of state: 



1 SF 
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= 2{A - DV^ 



(2) 
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The stability loss corresponds to an appearance of a 
non-trivial solution of the linearized equations of equilib- 
rium. Indeed, just at the brink of instability the system 
is in neutral equilibrium, if the non-linear terms are not 
taken into account. We shall look for this non-trivial 
solution for the x-component of the displacement vector 



u{y,z) 



(3) 



It has been mentioned above that the same result is ob- 
tained when the x and y axes are interchanged. 

Let us recall firstly that, while studying the equilib- 
rium of an elastic system, one has to consider separately 
the net strains, i.e. the homogeneous part of the strains, 
and their inhomogeneous part. The point is that the 
homogeneous part, which is described by six indepen- 
dent components of the (mean) strain tensor defines the 
change of the volume and the shape of the sample. The 
inhomogeneous part is better described by three indepen- 
dent components of displacement vector, since the use of 
the spatially varying strain components is usually less 
convenient as they should satisfy the additional (elastic 
compatibility) conditions. 

One can see immediately that in the present case there 
is no loss of stability with respect to homogeneous defor- 
mation. Indeed, if our sample (film plus substrate) has a 
homogeneous deformation (mJj,) , the energy of this de- 
formation would be 



F - {2AI + 2^L) «) 



(4) 



where L is the thickness of the substrate, and it would 
be infinite in the present case, L = oo. 

The inhomogeneous part of the strain should satisfy 
the equation of local equilibrium 



dy dz 
We shall apply the Fourier transform 

u{y, z) = / Ufe (z) exp {iky) dk 



(5) 



(6) 



and look for the appearance of the non-trivial solution of 
the equation of equilibrium (^ for a given wave vector k. 
We shall then determine the k where the stability sets in 
first (at highest temperature). Evidently, this would be 
the point where the symmetric phase looses its stability. 
For the ferroelastic (0 < z < Z) one obtains 



(Jxy [k, z) = 2AkUxy (k) = ikAkUk (z) , 



(7) 



where Ak = A + Dk"^, and 

^ duk 
o-xz (k, z) ^ A^-^- 

Equation (^ in the film {0 < z < I) takes the form 
(fuk Ak 



dz^ 



'^fc 7 2 n 

— k Uk = 0. 



(8) 
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For the substrate (— oo < z < 0) one has 



dz2 



k Uk — 0. 



(10) 



At the stress-free surface 

(z — I) the condition <Txz {k,l) — should be satisfied, 
or, taking into account Eq. (ph , 



duk 
dz 



\z=l= 



(11) 



At the interface (z = 0) the displacement Uk (z) and the 
stress axz (k, z) should be continuous, and the stress 
should vanish at z — > — oo. 

Let us first consider the case of < 0, which would 
correspond to a loss stability of the symmetric phase. 
The solution of Eqs. (|), (0) is 



Uk{z) = Fcosxkiz - I), 



(12) 



where — (k) /m, while for the substrate we have 



Ufe(z) = Gexp z. 



(13) 



The conditions at the interface give the equation for the 
existence of non-trivial solutions 



cotxfc^ = X- 



(14) 



This equation has a solution for > 0, while there is 
no solution for < 0, hence the loss of stability takes 
place for Ak < 0- For the region of interest we should 
obviously have x ^ 1, as we shall check later, and the 
approximate solution is 
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(15) 



(16) 



The minimum value of |A| = \^\c corresponds to the 
highest temperature where the instability sets in with 
k — km. '• 
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(17) 



1/2 

where usually (not on very "soft" substrates) {D/^) ' ~ 
dat with dat the interatomic distance (see, e.g., H). 
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Therefore, Xm ~ ^/kml ~ \/dat/l 1, and this con- 
firms our initial assumption that x is small. 

As a result, the loss of stability of the symmetric phase 
takes place at 



I 



(18) 



or, equivalently, at critical temperature, which is below 
that of a free sample: 



al 



(19) 



For displacive systems a fJ-/Tat, where Tat is the char- 
acteristic "atomic" temperature (see, e.g., [|ljl). One 
then finds that the lowering of the phase transition tem- 
perature is approximately 
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(20) 



Since usually Tat ^ (lO^ — 10"^) Tc^ one may expect a 
complete suppression of the second order transition in 
films with thickness of hundreds atomic layers, which is 
accompanied by the loss of stability of the symmetric 
phase. In order-disorder systems the shift of Tc is just 
(10~^ — 1)K. Certainly, one has to use the actual values of 
the coefficients D,a,fj, for a specific substance to decide 
if the second order transition is indeed suppressed in a 
given material. 

The stability gets lost with respect to inhoniogeneous 
strains u^y, u^z- The form of the "strain wave" in the 
film is 



Uxy = C sin {kraV + ^p) sin 
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Umz = C— — 7 COS [k„iy + (^) cos — , 
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(21) 
(22) 



where C is the constant, which remains undetermined 
within the linear theory. We see that the strain compo- 
nent u^y, which corresponds to the spontaneous strain 
in a free sample, is maximal near the free surface of the 
film, where it is easier to deform it. At the same time, 
the "accompanying" strain Uxz is concentrated near the 
film-substrate interface and is much smaller than the u^y 
strain for / ^ dat, because kml 3> 1. To find the ampli- 
tude C of the "strain wave" one has to take into account 
the non-linear terms in the equation of state. This is a 
difficult technical problem, and it is not addressed here. 
We shall just mention that at lower temperatures the si- 
nusoidal structure transforms into a "domain-like" struc- 
ture with the amplitude close to the value of the sponta- 
neous strain in a free crystal. One can interchange x and 
y to obtain an alternative variant of the "strain wave" . 

One can check that an inhomogeneous sinusoidal struc- 
ture appearing at a second order phase transition because 
of the boundary conditions is a precursor of a domain 



structure or, more precisely, it is the form of the do- 
main structure close to the phase transition, where the 
widths of the domain walls and the domains themselves 
become comparable (see, e.g. [|5j). Indeed, the domain 
wall width, Wc, is 



Wc. 
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^l/4;l/2 
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2\A\^ 2i/27ri/2^i/4 2fc„ 



(23) 



We see that indeed Vc is comparable with period of the 
"strain wave" at the temperature of the instability. 

We have found in Ref. that the usual values of the 
material coefficients involved the period a of the domain 
structure do not depend on temperature, and the domain 
width is given by the formula 



47r3 
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1/4 



(24) 



This is only about 50% more than the period of the 
"strain wave" appearing at the transition. In other 
words, the evolution of the "wave" into the domain struc- 
ture mainly consists of changing the shape of the wave 
(from sinusoidal to the piecewise domain-like) without a 
considerable change of its period. Since the period does 
not change much, the domain width becomes smaller 
than the domain wall width within the temperature inter- 
val about Tf. — Tcl- In a special case, when fj, is anoma- 
lously small and the film is thin enough, the period of 
the domain structure is much larger than that given by 
Eq. ( |2^ ) (see Ref. Q). At these conditions one should 
observe a more drastic evolution of the "strain wave" . 



III. STRAIN Uxx - Uzz AS THE ORDER 
PARAMETER 

We shall now consider the film where the soft mod- 
ulus corresponds to the strain Uxx — Uzz, Fig. 2. Once 
more, we assume that the elastic properties of the sys- 
tem are those of an isotropic medium with respect to all 
other strains. We shall see later that the character of the 
phase transition changes drastically if the misfit strain 
of the form assumed in Sec.l is nonzero. But we first 
consider the case where such a misfit is absent, to better 
appreciate its role. 



A. Domain pattern at the absence of extrinsic 
(permanent) misfit 

The same reasoning as in Sec. 2 leads to the conclu- 
sion that the soft wave vectors make 45° angle with the 
film plane. With the same reservations as in Sec. 2, the 
Landau free energy for the film substance is now: 
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FIG. 2. Schematic of the c/a/c/a "strain wave" and the 
domain structure (middle and bottom panels) in thin epitax- 
ial ferroelastic film with the spontaneous strain u'^^ — u1^. 
The symmetric phase looses its stability at some tempera- 
ture Tci < Tc, where Tc is the critical temperature in the 
free film, and forms the "strain wave" (middle panel). In the 
presence of the "extrinsic" misfit between the parent phase 
and the substrate (top panel) the "strain wave" evolves into 
the c/a/c/a domain structure with the period 2a — ai + 02 
at low temperatures in a film of thickness I. 
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where the critical strain Uc — Uxx — Uzz , A and /i are the 
Lame constants, and we shall consider the structures, 
which are homogeneous along the y axis. For the free 
energy density in the substrate we put yl = /i and omit 
the gradient term. 

Therefore, for the film we have 



cfxx = {X + H + A- DV^)uxx + (A - 
azz^ {>^ + ^^■ + A- DV^)uzz + (A - 

(yxz = '^^Uxz, 



- fi - A + D\'^)uz 
fj, ~ A + DS/'^)ux 



O-xx = (A + 2^i)Uxx + AWz2 
= (A + 2ll)Uzz + >^Uxx 



(30) 
(31) 



We can show that there is no loss of stabihty with respect 
to homogeneous strain. For inhomogeneous strains we 
have to satisfy the equations of equilibrium 



dx 
d<Jzz 



dz 



dz 
dx 



(32) 
(33) 



We have already mentioned in the previous section that 
the instability of the symmetrical phase is signaled by 
appearance of the non-trivial solutions of these equations 
(with the appropriate boundary conditions). 

We shall reduce the equations of equilibrium to a sys- 
tem of ordinary differential equations with the use of 
the Fourier transformation of the functions Ux {x, z) and 

Uz (X, z) 



(34) 



(z) [x, z)= I Ux(z) (fc, z) exp {ikx) dk 



We obtain the following equations for the Fourier com- 
ponents of the strain vector in the film (0 < z < ^) 



du. 



d^u^ 



fc^ (a + 1 - s^) Ux + ik{a + l + s^)'^ + H- ik 
(a -f 1 - 



duz 



dz 



dz ' dz^ ' dz 

+ik(a + l + ,s2) ^ - k^z + ik^ = 
^ dz '^■^ 



0, (; 



dz 



where a = A//i, — — {^A + Dk^^ / fx. The general solu- 
tion of this system reads 



^Cj-exp {xjz) , 




^(-l)^Qexp(x,z) 
(37) 



(38) 
(39) 



and we have kept only the first two terms of the expansion 
of Xj in terms of s ^ 1. For the substrate (— 00 < z < 0) 
the equations of equilibrium take the form 



du 



dz^ 



duz 



dz 



(27) -P{a + 2) Ux + ^ka-^ + -j-^ + ik-^ = 0, (40) 
(29) + 



-ika^^euz+ik^^O. (41) 
dz dz 



while for the substrate we obtain 



The solution that vanishes at z 
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(Di + D2Z) exp (kz) 



D1+D2 



(a + 1) fc 



exp {kz) 



(42) 
(43) 

Two conditions at 

the free surface {gzz — <^xz = at z = Z) and four con- 
ditions at the interface (continuity of Ux^ Uz, Czz, <^xz) 
provide six algebraic equations for Cj, Dj. This system 
of homogeneous equations has a nontrivial solution if the 
determinant of its coefficients is zero. After a straightfor- 
ward algebra, and again assuming that s 1, one finds 
the condition for nontrivial solution, 



cos 2kls 



1 



[a + iy (a + 2) 



(44) 



For < one has to replace the cos by cosh. It is clearly 
seen that there would be no solution for < 0, i.e. it is 
the case > that corresponds to the stability loss. 

Equation ( ^ ) refers to the stability loss for a given 
k. To find the actual stability loss one has, as in the 
previous section, to find the k vector that corresponds 
to the first instance where the instability sets in, i.e. at 
the highest temperature below Tc- To find it, one has to 
minimize the right hand side of the equation 



1^1 = Dk^ 



ii(a + l) o 
■ arccos 



4fc2/2 (q, + 2) 



(a + ly {a + 2) J 
(45) 



with the result 
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(a + iy{a + 2) 



We see that the condition {k„i) ^ dat/l ^ 1 is satisfied 
as long as / ^ dat- One sees that up to a numerical factor 
of order unity the result coincides with that given by Eq. 
(0). There is no essential changes if one puts A = 
(a = 0). We find the form of the "deformation wave" 
arising at the point of stability loss just for this case. 
Assuming once more that (km) «C 1 one finds: 



Uxx ix,z) = -Uzz (x^z) 



(48) 



= C cos [fc,„ (x ± z) + (p)] sin — , (49) 
Uxz = C —— sin [km {x±z) + Lp] cos — (50) 

where, once more, we obtain two possible "strain waves" 
corresponding to two possible periodic domain systems, 
from which only one can actually materialize. 

We found in Ref. that if the spontaneous shear 
strain appearing at the transition constitutes the only 



misfit strain existing in the system, then far from the 
transition (2 |j4| = ^) the period of the domain structure 
is given by the formula 
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(51) 



While obtaining this formula, the surface energy of the 
domain walls was estimated as D^/^Uq//^/^, where Uq is 
the spontaneous strain in a free sample. Such an estimate 
is reasonable for a proper ferroelastic. Comparing Eqs. 
( pT|) and (j^j we see that for the case in question the 
period of the domain structure far from the transition 
is close to its period at the transition where the domain 
structure is a sinusoidal "strain wave" . It is quite similar 
to what we have seen in Sec. 2. 



B. Effect of extrinsic misfit 

The properties of the domain pattern change substan- 
tially, if there is a misfit strain in the temperature range 
corresponding to the symmetric phase. This is usually 
the case, since the lattice constants of the film and the 
substrate are usually different. Therefore, the film is 
strained (with respect to the free film) at all tempera- 
tures. Since any misfit makes x- and z-axes inequivalent, 
the strain in the film includes the "symmetry breaking 
strain" , u^x—Uzz- In other words, the "symmetric phase" 
is no longer symmetric. Yet, a phase transition is possi- 
ble in this system. But, this is a transition between the 
single-domain and multi-domain states. It was studied 
[4:^ several authors |p|-[7|,p|JTl|] . It went unnoticed before 
that this transition is qualitatively different from the dis- 
/ ,X\ussed above: it is of the same type as the commensurate- 
^ incommensurate transitions and/or the transition be- 
tween the homogeneous and vortex states in type II su- 
perconductors. In other words, there is no loss of stabil- 
ity of the homogeneous (single-domain) phase at such a 
transition. 

The theory of this transition |ll| was developed for an 
extreme case which is opposite to the one considered in 
this paper. Indeed, the theory assumes that the width 
of the domain walls is negligibly small compared to all 
other length scales in the problem. We have found pl| ] 
that at the transition the period of the domain struc- 
ture diverges, while the width of the minority domains 
remains finite and equal to (27rZ)^^^ (£)//i)^^'' (in the no- 
tations of the present paper and neglecting a logarithmic 

factor). This length should be larger than the width of 

1/2 



the domain walls (Z?/ |yl|) ' , i.e. the inequality 



1^1 » 



(52) 
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should be obeyed. Comparing with Eq.(p^, we conclude [15] E.V. Chensky and V.V. Tarasenko, Sov. Phys. JETP 56, 
that the crossover between the two types of transitions 618 (1982) [Zh. Eksp. Teor. Fiz. 83, 1089 (1982)]. 

occurs in the region where the domain structure should 
be treated within a non-linear and non-local theory, but 
this is clearly beyond the scope of the present paper. 



IV. CONCLUSIONS 



We have considered two examples of proper ferroelas- 
tic transitions in thin epitaxial films. In the first one the 
ferroelastic strain Uxy breaks the symmetry in the film 
plane, which was not broken by the "extrinsic" misfit (ir- 
relevant to the that phase transition). There a "usual" 
(Landau-like) second order phase transition takes place 
with a loss of stability of the symmetric phase and forma- 
tion of sinusoidal "strain wave" . The period of the wave 
is close to the period of the domain structure, which will 
evolve at lower temperatures, except for a special case of 
very thin films and anomalously small non-critical elastic 
modulus fj,. In the second case of the ferroelastic Uxx—Uzz 
strain the results are similar when the "extrinsic" misfit 
(mismatch between the lattice parameters of the sym- 
metric phase of the film and the substrate) is absent. If 
the "extrinsic" misfit is present, the phase transition be- 
comes the (topological) transition between a single- and 
a multi-domain states. The crossover between the tran- 
sitions of the two types remains an open issue. 
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